On the role of electroweak bremsstrahlung for indirect dark matter signatures 
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Interpretations of indirect searches for dark matter (DM) require theoretical predictions for the 
annihilation or decay rates of DM into stable particles of the standard model. These predictions 
include usually only final states accessible as lowest order tree-level processes, with electromagnetic 
bremsstrahlung and the loop-suppressed two gamma-ray line as exceptions. We show that this 
restriction may lead to severely biased results for DM tailored to produce only leptons in final states 
and with mass in the TeV range. For such models, unavoidable electroweak bremsstrahlung of Z and 
W-bosons has a significant influence both on the branching ratio and the spectral shape of the final 
state particles. We work out the consequences for two situations: Firstly, the idealized case where 
DM annihilates at tree level with 100% branching ratio into neutrinos. For a given cross section, 
this leads eventually to "minimal yields" of photons, electrons, positrons and antiprotons. Secondly, 
the case where the only allowed two-body final states are electrons. The latter case is typical of 
models aimed at fitting cosmic ray e~ and e + data. We find that the multimessenger signatures of 
such models can be significantly modified with respect to results presented in the literature. 



PACS numbers: 95.35.+d, 95.85.Pw, 98.70.Sa 

I. INTRODUCTION 

Despite the numerous cosmological and astrophysical 
indications for the presence of nonbaryonic dark matter 
(DM), the particle nature of DM has yet to be identified. 
One strategy towards DM "detection" is to search for 
its self-annihilation (or decay) products in our Galaxy, 
provided that the annihilation cross section (av) (or the 
decay rate) is large enough. Restricting the space of 
candidates to weakly interacting massive thermal relics, 
(av) (at freeze-out) is fixed by the DM abundance while 
the DM mass mx lies in this class of models typically 
within one or two orders of magnitude off the weak scale 
mw- Nevertheless, considerable model-dependence re- 
mains because of our ignorance of the final states pro- 
duced in the annihilation process, which vary from model 
to model. 

Under the sole hypothesis that massive dark matter 
annihilates into standard model particles, a few years 
ago an interesting conservative upper bound on (av) was 
derived using data on the least detectable final states, 
namely neutrinos [1, 2]. (Similar considerations apply to 
decaying particles [3] , although we will not mention this 
further.) It is natural to ask if these conservative bounds 
can be improved or consolidated further. It was shown 
previously in Ref. [4] that electroweak bremsstrahlung 
leads to a break-down of perturbation theory and a non- 
negligible branching ratio in electromagnetic channels for 
nix 3> "iff, even if at tree level DM couples only neu- 
trinos. This argument was used then in Refs. [5, 6] to 
derive constraints on the DM annihilation cross sections 
from diffuse gamma-ray observations, that turned out to 
be similarly restrictive as the original one from Ref. [1, 2]. 

More recently, the PAMELA collaboration published 
data showing an "anomalous growth" of the cosmic 



ray positron fraction [7], while the measured antipro- 
ton fraction agrees with expectations from simples mod- 
els [8] . Additionally, new data on the overall electron plus 
positron spectrum were presented, most notably from the 
Fermi space telescope [9, 10]. These data have prompted 
a plethora of models trying to explain the data by engi- 
neering relatively heavy (mx > 1 TeV) DM candidates 
with large annihilation cross sections (or decay rates) and 
small or vanishing branching ratios (br's) into hadronic 
final states. In the analysis of these models, the role of 
radiated W, Z bosons and their phenomenological impact 
has been generally ignored. 

The purpose of this article is to revisit the topic of 
electroweak bremsstrahlung effects with several goals in 
mind. First of all, albeit the qualitative conclusions of 
Refs. [5] and [6] agree, their results differ quantitatively 
Here, we repeat this calculations analytically and cross- 
check them against numerical results using Madgraph. 
Second and more important for its recent phenomeno- 
logical appeal, similar effects also arise when the tree- 
level final state is a charged lepton £ + £~, and therefore 
we present an analogous calculation for an e + e _ pair as 
final state. Third, besides photons, the fragmentation 
of the emitted W and Z bosons produces also electrons 
and positrons, protons and antiprotons, as well as neu- 
trinos and antincutrinos (hereafter, simply "neutrinos"). 
Hence we can use observations of different cosmic species 
to derive complementary constraints on DM annihila- 
tions. Finally, we include in the present analysis new 
data on the antiproton [8] and positron fraction [7] from 
the PAMELA satellite as well as on diffuse gamma-rays 
from the Fermi Gamma-ray Space Telescope [11]. 

This article is organized as follows: In Sec. II 
we review first some general considerations about the 
favoured DM annihilation modes and then our analyt- 
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ical calculation of the branching ratio for electroweak 
bremsstrahlung. In Sec. Ill we describe the spectra of 
secondaries found numerically, while Sec. IV is devoted 
to the constraints from present observations. We con- 
clude in Sec. V. 



II. ANALYTICAL CALCULATION OF B.R.'S 
INTO W,Z. 

A. General considerations 

It is interesting to ask oneself under which circum- 
stances it is possible to produce unsuppressed or even 
dominant final states consisting of (possibly light or 
massless) leptons. In the following, for the sake of self- 
consistency, we elaborate on some of the considerations 
reported in the Appendix of [6] and in [12]. 

On general grounds, the L-th partial wave contribu- 
tion to the annihilation rate of two heavy, non-rclativistic 
particles moving with relative velocity v is suppressed 
as v 2L . The virial velocity in our Galactic halo is only 
v ~ 10~ 3 . Typically, only the L = partial wave results 
in observable rates for indirect DM detection today. Ad- 
mitting L = 1 final states, terms proportional to v 2 and 
terms of similar magnitude that are chirally suppressed 
as oc (rrif /mx) 2 , where vtif denote the mass of fcrmions, 
enter the corresponding annihilation rate: This is no- 
tably the case for annihilations via an axial vector cur- 
rent, J PC = 1 ++ . (We use here the spin-parity notation 
J PC , where J is the total — orbital plus internal — spin 
quantum number, P the parity and C the charge-parity). 

Such a scenario has two classes of phcnomcnological 
problems: First, any observable effect now should pro- 
duce huge effects at early times, when v was larger. Sec- 
ond, multi-body final states, e.g. with e.m. or weak bo- 
son radiation emitted from one fermion, proceeds unsup- 
pressed as shown for the case of the axial vector mode in 
Ref. [13]. In this case, bounds from photon and cosmic 
ray antimatter searches typically rule-out these models. 
This provides a general argument to focus on L = 0- 
modes only, as we do in the following. 

For a spin- 1/2 Majorana fermion candidate, the par- 
ity is P = (— 1) L+1 and the charge parity state C = 

(-l)o r bital(-a S ptnVl)a„ticom m . = (-l) i+S , where the 

spin S = 0, 1 for a fermion pair. Since L = is the 
only acceptable choice, it follows that P = — 1. Also, the 
Majorana nature of the particles implies that it must be 
even under C, hence S = 0. The only state selected is 
thus the pseudoscalar h . 

Another class of popular DM candidates are scalars. 
Again, since we are in the non-relativistic limit, the only 
unsuppressed annihilation state is the scalar one, with 
L = — > P = — 1. Either a fundamental DM scalar or 
pseudoscalar would lead to a scalar singlet. The state is 
clearly even under C-symmetry. Thus, J PC = ++ , i.e. 
only a scalar is allowed. A spin-0 two body particle states 
is also allowed in case of Dirac and vector particles as 



DM candidates, although it is then not the unique choice 
leading to L = annihilations [12]. 

In summary, for the following calculations we can 
adopt an effective field theory approach and define 
the initial state as T> coupling with SM neutrinos v 
as — \\T>Vv or the analogous pseudoscalar coupling 
— iA W^v. Basically without loss of generality, this is a 
viable way to obtain on a phenomenological level a tree- 
level coupling to neutrinos only. Note that the exact 
flavour-structure of the coupling is irrelevant, since the 
two large neutrino mixing angles lead after oscillations 
to a 1 : 1 : 1 mixture of flavours. A similar coupling will 
be adopted for the electron final state case as well. 

B. Explicit calculation 

As discussed above, we can replace e.g. the annihilation 
process XX — > VvZ with the decay T> — > VvZ choosing 
rriT> = 2m x- In particular, this replacement becomes 
exact, if the annihilation process is mediated by a scalar 
particle, 

_ a (XX -> VvZ) _ r(D -> VvZ) 
Z ~ a(XX -> Vv) ~ T(V -> Vv) ' ( ' 

To be specific, we choose a scalar coupling Ci = — XDVv 
between neutrinos and T>. The tree-level decay rate of 
T> in one massless neutrino flavour is then T(T> — > vv) = 
A 2 mxi/(87r). Note that this is a factor two smaller than 
the decay rate used in Ref. [6]. Also note that in the 
above expression v is implicitly assumed to be a Dirac 
field, hence equal populations of active left-handed neu- 
trinos vl and "sterile" right-handed states vr are pro- 
duced. While there exists no evidence for the existence 
of vrs up to now, this choice leads to the most conser- 
vative bounds: Since the "sterile" right-handed neutrino 
states vr do not participate in electroweak interactions, 
the br's of the radiative channels would be enhanced re- 
stricting the tree-level coupling to vl only. 

The two diagrams describing the electroweak 
Bremsstrahlung process and our notation for the 
momenta are shown in Fig. 1. The corresponding 
Feynman amplitudes are 

M 1 = u( Pl ) -^-rP L -J- e„(p 3 ) (-iA) v( P2 ) , (2) 
2c w <fr l 

M 2 = u( Pl ) (-iA) i- -±VPl e„(ps) v(p 2 ) , (3) 

where we have neglected neutrino masses, m„ = 0, as 
always in the following. Furthermore, g is the weak cou- 
pling constant, c\y = cosflw and Pl/r = (1 ± 7' 5 )/2 
are the left/right-chiral projection operators. Using as 
short-hand notation for the polarization tensor of a mas- 
sive gauge boson P M „ = -g^ v + p^Pn, / f m% as well as 
Ki ~ [(9^)/(2cwqf)] 2 , we obtain summing over spins 

JAjA\ = K.P^Tt^Pl^^Y) = KyP^B^ . 

(4) 
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The interference term vanishes for m v = and thus we 
can set \Ai\ 2 = 2\A4i\ 2 in the calculation of the decay 
rate. Evaluating the trace in the rest-frame of T> gives 
with p = (mp, 0) 



P^Bf _ 2 {ml + m 2 z - m{ 2 ) 



m 



v 



m\ 



(5) 
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{m% + m| - mf 2 - m§ 3 ) --z 

where we introduced m?^ = (p.; + Pj) 2 = (p — Pk) 2 [14]. 

The differential decay rate for a general 1 — > 3 decay 
process, 

1 1 



dr = 



(2tt) 3 32m| 



\M\ dmj 2 dm 



2 

23 ' 



(6) 



can be integrated exactly using the limits [15] 
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The resulting ratio Rz is 
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with yz = p 2 /m 2 z = Am 2 x /m 2 z . The ratio Rw follows 
from Rw(yw) = ^c^Rziyw)- Our ratios have the same 
dependence on yz,w as those found by the authors of 
Ref. [6], but are a factor four larger. A factor of two is 
explained by the difference in the tree-level decay width, 
while the other, given the agreement of our analytical 
results with numerical ones presented in the following, is 
attributed to a missing algebraic factor in the calculation 
of [6]. 

For the electron final state case, the ratio Rw writes 
similarly as above, but the different structure of the cou- 
pling with Z (involving also Pr) does not lead to a sim- 
ilar simple expression for Rz', in particular, interference 
terms do not vanish in the unitary gauge. The latter 
contribution has been thus calculated only numerically. 

In our calculation of Rz we have included only fi- 
nal state radiation (FSR) neglecting virtual state radi- 
ation (VSR) "from the decay vertex" as well as initial 



state radiation (ISR). The latter two can be only calcu- 
lated within specific models or assuming a separation of 
scales such that an effective theory approach can be used. 
The separation between three classes of bremsstrahlung 
is gauge-dependent and thus strictly speaking meaning- 
less. However, in cases where one of the three radiation 
mechanisms dominates and no cancellations are present, 
this separation is useful. For the electromagnetic case, 
the relative importance of VSR and FSR is discussed 
e.g. in Ref. [13, 16, 17]. In some cases, VSR is impor- 
tant since, depending on the exact spin structure of the 
interaction and the particles involved, the helicity sup- 
pression present for certain two body final states can dis- 
appear. This is the opposite limit to the one we are dis- 
cussing here, since we are working within the ansatz of 
dominant neutrino or electron two-body final states. For 
our purpose, it is reasonable to conclude that neglecting 
the (model-dependent) ISR and VSR provides at most a 
slightly conservative evaluation of the visible final state 
channels. 

Also, it has been shown in [18] that, in addition to the 
Z and IV-strahlung processes, loop-induced branching 
ratios into charged leptons (and, we note, more generally 
into charged fermions including quarks) arc unavoidable 
in these scenarios. Similarly to the VSR, these effects de- 
pend on the UV completion of the model and we do not 
discuss this further. Typically, one expects these branch- 
ing ratios to be at the percent level, and thus mainly 
important for relatively light DM particles, mx < mw- 
In any case, loop processes as additional source of pho- 
tons and anti-matter would strengthen the limits derived 
in the following. 

Let us discuss now the behavior of the bremsstrahlung 
amplitudes in the limit m-p 3> mz- The similar cal- 
culation of the three-jet rate in QCD with the emis- 
sion of a massive gluon leads to -Rqcd oc In 2 (m^/m 2 ), 
and in this process the longitudinal gluon does not con- 
tribute as it is coupled to a conserved current. By con- 
trast, the axial vector current is not conserved and it is 
indeed the longitudinal polarization of the clcctroweak 
gauge bosons that is responsible for the quadratic mass 
dependence, Rt oc mf, [19]. This difference is also the 
main reason for the discrepancy between the results of 
Ref. [5] and [6]. In the former work, the probabil- 
ity for electroweak bremsstrahlung was calculated using 
the Feynman-'t Hooft gauge and as tree level process 
v'v' — > vv. Hence, additionally to the differences ex- 
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pected from using a different current and coupling struc- 
ture mediating the annihilations, the slower rise of the 
br's (which also holds in SUSY, see below) found in [5] 
can be attributed to the absence of the longitudinal po- 
larization of the gauge bosons. 

Unitarity requires e.g. that the annihilation cross sec- 
tion behaves as (<rv) oc m^ 2 . While this bound is re- 
spected by the tree-level cross section, the Rf oc mf, 
dependence leads to a violation of pcrturbative unitar- 
ity if an additional W or Z is emitted. The decoupling 
property of supersymmetry [20] guaranties that such a 
quadratic term is absent, even in the presence of soft 
masses 3> myy- An explicit calculation in the MSSM 
show that the ratio of e.g. cr{xx ~ * e + e~ Z)/a(xx ~ * 
e+er) ~ 0.03 for a lOTcV neutralino. Hence the MSSM 
is an example for a scheme where initial, virtual and fi- 
nal state radiation from all sub-processes in the s, t and 
u channel is arranged in such a way that their leading 
contributions cancel. On the other hand, in such a the- 
ory quarks and gauge bosons are produced at tree-level 
in two-body final states, so that the yields of secondaries 
other than leptons are typically large. 

Similarly, the above considerations can be circum- 
vented if the DM does not decay directly into SM par- 
ticles, rather trough some light state Y , so that W, Z in 
the final state are kinematically inaccessible. These mod- 
els however introduce other light, metastable degrees of 
freedom, which violate from the beginning our "effec- 
tive theory" approach. On the other hand, for a chain 
XX — > nY — > SM, provided that my > mz, one might 
still apply the above arguments with the new ratios writ- 
ten now as R' z w = nR(m Y /m 2 z w ), although the distri- 
butions of secondaries will be different, of course. 



III. NUMERICAL RESULTS 



We have calculated the total cross section XX — > 
v e e W T also numerically, using Madgraph [21] and al- 
lowing electroweak bremsstrahlung only as FSR. More 
precisely, we used the MSSM model, choosing X as the 
lightest neutralino and the scalar higgs h as intermedi- 
ate particle as well as switching off all diagrams except 
the one corresponding to FSR. We denote the ratio of 
the bremsstrahlung and the tree-level processes again as 
Rw, 



Rw = — 

02 



a(XX -» v e e+W~) + a(XX 



-w+) 



a(XX — ► v e v e ) 



(10) 

In Fig. 2, we compare our analytical result for Rw with 
the numerical results obtained with Madgraph, finding 
excellent agreement. 

Note that the perturbative results for the 
bremsstrahlung cross sections become unreliable al- 
ready at mx Si 1 TeV. Above this energy, processes with 
n > 3 (treating W, Z as stable) particles in the final 
state become important and an electroweak cascade 
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FIG. 2: Our analytical result for the ratio Rw (solid line), 
compared to Rw from ref. [6] (dashed line) and the numerical 
results from Madgraph (points). 



develops [4]. On the other hand, each individual sub- 
process is suppressed by a Sudakov factor compared to 
the perturbative result, avoiding a "blowing-up" of the 
total annihilation cross section. 

In the next step, we feed the events generated by Mad- 
graph into HERWIG++ [22] and generate the energy 
spectra of e ± , vi,Vi, p, and p produced as secondaries in 
W and Z decays. The obtained energy spectrum dN/dx 
are normalized to the total cross section, crtot = (J 2 + c.3 , 



dNi _ dcr(XX -» i + all) 
da; 



Ctot 



d.r 



(11) 



where /j is the multiplicity to produce particles of type 
i with energy E — xmx- Therefore the evolution of 
dN / dx with energy becomes much slower for mx ~ few 
TeV, when cr tot « 03 » a%. Since the shape of the 
fragmentation function dNi /dx changes only logarithmi- 
cally, the omission of 2 — ► n > 3 processes has only a 
minor impact on dNi/dx. Moreover, the chosen normal- 
ization gives the correct number of secondaries Ni per 
annihilation also for 03 3> 02. The resulting fragmen- 
tation functions dN/dx as function of x are shown for 
m x = 300 GeV and m x = 3 TeV in Fig. 3. Note that i) 
there is an additional contribution 5(1— x) 02/otot wn i cn 
is not reported for clarity and ii) only half of the neutri- 
nos emitted at tree-level are "active" ones, vl. 

In the same way as described above for tree-level an- 
nihilations into neutrinos, we have calculated also the 
total cross sections of XX — > e + e~ and XX — > e + e~ Z, 
XX -> v^W^. Additionally, we added analytically 
photons from external electromagnetic bremsstrahlung, 



diV 7 
dx 



a 



Ps(x) In 



8(1 - X) 



(12) 
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FIG. 3: The fragmentation spectra dN/dx of electrons (long-dashed maroon), photons (solid red), neutrinos (dot-dot-dashed 
green), and protons (dotted black) in "only neutrinos" tree level annihilations for mx = 300 GeV (left) and mx = 3TeV 
(right). The monochromatic neutrino "spike" at x — 1 is not shown. Fermion labels indicate the sum of matter and antimatter. 



with the usual fcrmion-fermion-boson splitting function 



The resulting fragmentation functions dNi/dx are 
shown for mx = 300 GeV and mx = 3TeV in Fig. 4. 
Note that the relative importance of protons is largest for 
x ~ 0.1 where it is comparable with the secondary elec- 
tron flux. Also, brcmsstrahlung photons (absent for the 
neutrino case) provide the dominant yield only for sub- 
TeV DM masses, while for heavy DM particles, apart 
for the region x > 0.1 and the very small x, the pho- 
tons from W, Z fragmentation dominate. As above, there 
is an additional contribution to the e~ + e + case of 
2 8(1 — x) 02I o to t which is not reported for clarity, but of 
course included in the constraints derived below. Note 
also that, in contrast to the neutrino case, both helicity 
states emitted at tree-level contribute to the observable 
electron spectrum. 

IV. CONSTRAINTS 

In this section we provide a first application of the cal- 
culations reported above, deriving bounds on (av) vs. 
mx from measurements of the diffuse gamma ray flux, 
the antiproton fraction, the electron plus positron flux, 
the positron fraction and the limits on Galactic neutrino 
fluxes. We shall adopt a simple prescription, namely we 
consider a model as excluded when the DM contribution 
to the signal alone exceeds the maximal flux allowed by 
the data at one sigma. This is likely over-conservative 
since it is known that other more mundane astrophys- 
ical sources contribute to or even dominate the fluxes. 



For easing the comparison with other papers, we adopt 
for the smooth DM mass density p sm a fiducial Navarro- 
Frcnk- White profile [23] 



where a = 20 kpc is the characteristic scale below which 
the profile scales as r -1 . Following the new determina- 
tion in [24], we choose p Q = 0.39GcV/cm 3 as the DM 
density at the solar distance from the (Galactic center) 
GC; the latter is estimated to be r Q sa 8.33 kpc [25]. 
Since we shall refer to high-latitude fluxes (for the photon 
case) or diffuse signals (for the charged particles cases) 
the exact DM profile towards the GC is not essential. 
For neutrinos, we shall limit this dependence adopting a 
quite large cone size (see below). 

Note that additional bounds can be derived by focusing 
on more specific location (as the GC) and/or by looking 
at other channels produced as secondaries of leptonic en- 
ergy losses (inverse Compton, synchrotron radiation,. . . ). 
It is beyond the scope of this paper to provide an ex- 
haustive account of indirect bounds. Here we only note 
that most of them depend more on the properties of the 
Galactic medium, and that further constraints can only 
strengthen the results presented here. 
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FIG. 4: The fragmentation spectra dN/dx of electrons (long-dashed maroon), photons (solid red), neutrinos (dot-dot-dashed 
green), and protons (dotted black) in "only electron" annihilations for mx = 300 GeV (left) and mx = 3TeV (right). Also 
shown is the gamma-ray yield from bremsstrahlung only (short-dashed red). The monochromatic electron "spike" at x = 1 is 
not shown. Fermion labels indicate the sum of matter and antimatter. 



A. Gamma-rays 



similarly to Eq. (15) as 



The differential flux of photons from (self-conjugated) 
dark matter annihilations is 



dv 7 Mp| 

dE 8T:m x 



ds 



n.o.s. \ Po 

(15) 

where r(s,ip) = (r| + s 2 — 2r Q scosip) 1 ^ 2 , ip is the an- 
gle between the direction in the sky and the GC and s 
the distance from the Sun along the line-of-sight (l.o.s.). 
In terms of the galactic latitude b and longitude I, one 
has cos ip = cos b cos I . Particle physics enters via the 
DM mass mx, the annihilation cross section (av), and 
the photon differential energy spectrum dN 1 /dE per an- 
nihilation. Since here we are focusing on the galac- 
tic diffuse emission rather than that from the GC, the 
residual uncertainties which arc introduced through the 
choice of the DM profile are negligible for our discus- 
sion. We shall compare the calculated flux with the 
high-latitude residual "isotropic" emission whose prelim- 
inary data have been presented by the Fermi-LAT col- 
laboration [11], whose upper limit at ~ 1 a in the range 
0.1-50 GeV can be roughly represented by E 2 dN/dE 
1.5 x 10- 2 (£/0.1GeV)-°- 45 ' 



'GeVm- 2 s _1 sr- 



B. Antiprotons 

Accounting for diffusion, the flux of antiprotons at 
Earth is isotropic to a high accuracy. By neglecting en- 
ergy losses and the so-called "tertiary" component (which 
is mostly relevant at low-energies) the flux can be written 



d^^ 
dE 8 7r mi 



(16) 



where dNp/dE is now the antiproton injection spectrum 
per annihilation, while the integral along the line of sight 
is effectively replaced by a function Fp{E) which encodes 
the dependence of the flux from astrophysical parame- 
ters (see [26] for a derivation of the above formula and 
an explicit expression of Tp ) . For the present illustrative 
purposes, we adopt the fit of this function calculated nu- 
merically in [27], for the reference NFW model and the 
"intermediate" choice of propagation parameters (see [27] 
for details). 

In order to compare with the p/p ratio provided in [8], 
we normalize the above calculated flux to the observed 
proton flux which we take from [28] to be 



— ? — i 
m s sr 



$; bs (£) = 1.4 x 10 4 (£/GeV)-'' 7 GeV 

(17) 

where we accounted for a proton contribution of 79% in 
the overall cosmic ray flux at the energies of interest. 



Electrons and Positrons 



Compared to the above case, the main difference con- 



cerning the calculation of the flux of e = (e 



") at the 



Earth starting from the injected parameters is due to the 
relevance of energy losses. A semi-analytical form can be 
derived [29, 30], which yields 



dE 8- 



(18) 
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with 



1 



^k(E)b(E) J E 



d£^(£)h e (£). (19) 
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In the above equation, b(E) « 10~ 16 (E/GeV) 2 GeVs" 1 
encodes energy losses and the function h e (E) depends 
on halo parameters as well as propagation ones. As for 
the case of antiprotons, we adopt the fit for this function 
provided in [27] for the NFW case and intermediate prop- 
agation parameters. Note that the positron flux amounts 
to half of the flux reported above. 

From the overall electron plus positron spectrum we 
know from the Fermi measurement that [9, 10] 



<f>° hs (E) = $ 



obs 
+ +c- 



(E) 



(20) 



(175.40 ±6.09) 



E 
GeV 



-3.045±0.008 



1 



GeV m 2 s sr 



represents a good fit of the data between ~ 20 GeV and 
1 TeV. We shall then require that: i) $ c /(2$° bs ) < 
f"+(E), where f"+(E) is the upper value of the data- 
points presented by PAMELA at energies between ~ 20 
and ~ 100 GeV. ii) That $ c < $° bs > u P, where $° bs < u P is 
the flux obtained for the upper value of the flux normal- 
ization and the hardest spectrum reported in Eq. (20). 



D. Neutrinos 

For neutrinos, we use the Supcr-Kamiokandc 
bound [31] requiring that the induced flux of muon 
tracks from a cone of half-width 30° around the GC 
is below 1.6xl0~ 14 /s (we specifically use the slightly 
conservative prescription to use only events induced by 
neutrinos above 10 GeV energy). Also, we assume a 
1:1:1 mixture of flavours, which holds within factors 
of order unity independently of the production flavour 
ratios due to mixing. Finally, note that the width 
of the cone around the Galactic Center is such that 
the dependence of the flux on the exact DM profile is 
marginal for our purposes (within a factor 2 at most, 
see [2]). 



E. Results 

In Fig. 5 we summarize our bounds. Not surpris- 
ingly, for the case of an electron final state at tree level 
(top panel), the bounds from electron-related observables 
dominate at low energies. In particular, around mx ~ 
100 GeV, the positron fraction is saturated by a cross 
section less than one order of magnitude larger than the 
fiducial value for a thermal relic, (av) ~ 3 x 10~ 26 cm 3 /s. 
Equivalently, we expect a DM contribution above the 
10% level in some bins of the PAMELA data, even in 
absence of astrophysical boost factors due to a clumpy 
DM distribution. Note also that for the same mass range 
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FIG. 5: In the top (bottom) panel, we report the exclusion 
plot in the (av) vs. mx plane for the case of electron (neu- 
trino) final states at tree level. Regions above the different 
curves are excluded by the labelled CR observables as de- 
scribed in the text. 



the diffuse gamma-ray bounds fall, despite the a/ir sup- 
pression, within two orders of magnitude of the fiducial 
value, confirming the promising role of this channel for 
detection in case of "ordinary" final state br's and of 
more targeted searches (see e.g. [32] and refs. therein). 
Above mx ~ TeV, the limits derived from different chan- 
nels have roughly the same strength. Hence saturating 
the electron flux by DM annihilations will lead to ten- 
sions with the antiproton bound, and eventually with the 
gamma-ray and neutrino bounds too. 

For the neutrino final state the best limit comes from 
neutrino observations, but for the region of a few TeV 
where the antipronton fraction provides a better con- 
straint. Also, in the same range the other channels lead 
to constraints looser by less than one order of magnitude, 
with the gamma-ray being the most competitive one. 
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In deriving the above bounds, we assumed conserva- 
tively that the products of the electron energy losses 
are not observed otherwise, while it is known that they 
can lead to potentially interesting signatures in "softer" 
gamma-rays via inverse Compton scattering, for exam- 
ple. Also, we ignored gamma-ray observations closer to 
the GC, while it is estimated that gamma-ray constraints 
at intermediate galactic latitudes or towards the GC are 
stronger than the ones presented here (see e.g. [33-35].) 
Finally, enforcing the HESS constraints [38] on the com- 
bined e~ + e + + 7 flux at energies E « 0.6-4 TeV might 
produce — depending on the propagation parameters — 
bounds a factor of a few more stringent than the ones 
shown. Even with these caveats in mind, however, it it 
interesting to note that multi-TeV scale particles annihi- 
lating into neutrinos at tree-level might turn out to be 
similarly or even more constrained than particles annihi- 
lating into electrons. 

The comparison with the electron case shows that even 
for the largest masses the tree-level, monochromatic neu- 
trinos provide the leading constraint. This is essentially 
due to two facts: i) The cross section for muon produc- 
tion and the muon range continue to increase also for 
s 3> m 2 ^, partially compensating the m^- 2 suppression 
in the differential flux, ii) The lack of energy losses for 
neutrinos, which is instead a key factor for very energetic 
electrons. 

On general grounds, it appears safe to conclude that if 
a DM candidates in the few TeV mass range contributes 
a significant fraction of the CR electron/positron flux, 
then a non-negligible antiproton flux is produced as well 
that is close to the current bound. Especially in the light 
of the more refined measurements expected in the next 
years with AMS-02 [36] , this is an important signature to 
keep in mind. Further, for the same class of candidates 
one expects measurable signatures in both the gamma- 
ray and the neutrino flux, with interesting observational 
perspectives from the direction of the inner Galaxy. 

V. CONCLUSIONS 

We discussed the role of electroweak bremsstrahlung 
for indirect dark matter signatures, which has been typ- 
ically neglected in phcnomcnological studies heretofore. 
Our approach was to calculate the branching ratios and 
the fragmentation spectra diV^/dx of secondaries con- 
sidering electroweak radiation only from the final state. 



Therefore, our results are applicable mainly to models 
tailored to produce only leptons as final state and with 
DM mass in the TeV range. In particular, our analy- 
sis applies to several models trying to match features in 
e + e~ CR data with TeV-scale DM, but e.g. not to the 
benchmark case of neutralino annihilations in the MSSM 
(where final states containing hadrons or gauge bosons 
are anyway allowed at tree- level). 

An important phenomenological consequence of elec- 
troweak bremsstrahlung is its impact on the predicted 
photon and electron/positron spectra. Secondaries from 
W and Z decays can dominate in a certain x range the 
spectra, leading to changes in the normalization and the 
shape of the secondary spectra. Thus it is mandatory 
to assess the importance of these "corrections" in a spe- 
cific model under consideration, before one attempts to 
fit cosmic ray or photon data. Models viable at tree level 
may be ruled out by the more realistic spectrum or, vice 
versa, large boost factors required might be moderately 
loosened. 

On more general grounds, at least in cases where 
no new light particles are introduced in the spectrum, 
it appears that DM candidates in the few TeV mass 
range contributing a significant fraction of the CR elec- 
tron/positron flux should lead also to a non- negligible 
antiproton flux, close to the current bound. Addition- 
ally, for the same class of candidates we expect measur- 
able signatures in both the gamma-ray and the neutrino 
flux. These "multimessenger" signatures are not surpris- 
ing and have often been implicitly assumed in past phe- 
nomenological works on indirect DM detection. How- 
ever, it is interesting that the conclusion holds virtually 
unchanged also for heavy DM candidates engineered to 
produce only leptons as tree level final states. 
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